Slow Relaxation in a Constrained Ising Spin Chain: a Toy Model for Granular 

Compaction 
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We present detailed analytical studies on the zero temperature coarsening dynamics in an Ising 
spin chain in presence of a dynamically induced field that favors locally the ' — ' phase compared to 
the '+' phase. We show that the presence of such a local kinetic bias drives the system into a late 
time state with average magnetization m = — 1. However the magnetization relaxes into this final 
value extremely slowly in an inverse logarithmic fashion. We further map this spin model exactly 
onto a simple lattice model of granular compaction that includes the minimal microscopic moves 
needed for compaction. This toy model then predicts analytically an inverse logarithmic law for 
the growth of density of granular particles, as seen in recent experiments and thereby provides a 
new mechanism for the inverse logarithmic relaxation. Our analysis utilizes an independent interval 
approximation for the particle and the hole clusters and is argued to be exact at late times (supported 
also by numerical simulations). 



PACS numbers: 05.40.-a, 82.20.Mj 
I. INTRODUCTION 



Slow relaxation dynamics naturally occurs in systems 
with quenched disorder such as spin glasses and has re- 
mained a subject of long standing interest However, 
systems without quenched disorder such as structural 
glasses also exhibit slow dynamics. It is believed that 
the slow relaxation in the latter systems is due to ki- 
netic disorder, induced by the dynamics itself An- 
other important class of systems without quenched dis- 
order is granular material, where once again kinetic dis- 
orders are responsible for slow relaxation. In a recent ex- 
periment ||], a cylinder packed loosely with glass beads 
was tapped mechanically and it was found that the sys- 
tem gets more and more compact with time. However, 
the density p(t) compactificd rather slowly with time as, 
p(oo) — p(t) ~ l/log(£). How robust is this inverse loga- 
rithmic relaxation? Is it only specific to granular systems 
or does this also occur in other out of equilibrium systems 
in presence of kinetic disorders? 

In this paper we study, in detail, the effect of kinetic 
disorders on the relaxation dynamics of an Ising system 
quenched from a high temperature disordered phase into 
the low temperature ordered phase. In absence of ki- 
netic disorders, the dynamics of such a system is well 
understood [Q. As time progresses, domains of equi- 
librium low-temperature ordered phases (consisting pre- 
dominantly of up and down spins respectively) form and 
grow. The average linear size of a domain grows with 
time as l(t) ~ t = for zero temperature nonconserved dissi- 
pative dynamics. How does disorder, quenched or kinetic, 
affect this simple dynamics? The effect of quenched fer- 
romagnetic disorder on this phase ordering kinetics has 
been studied extensively Q. Essentially the quenched 
disorder tends to pin the domain walls leading to com- 
plete freezing at T = 0. However at small nonzero tem- 



peratures, the domains still coarsen via activated dynam- 
ics but extremely slowly as, l(t) ~ (logt) 1 / 4 The 
purpose of this paper is to explore the effect of kinetic 
disorder on the phase ordering dynamics. A shorter ver- 
sion of this paper with the main analytical results and 
their numerical confirmation has appeared elsewhere [EJ . 
Here we explore the dynamics in more detail and also 
elaborate its connection to other systems such as granu- 
lar material and reaction-diffusion systems. 

In this paper, we study the zero temperature dynamics 
of an Ising chain in presence of a specific type of kinetic 
disorder namely a dynamically induced magnetic field. If 
a small uniform external field (say in the down direction) 
is put on in an Ising system following a rapid quench from 
infinite temperature to T = 0, the system rapidly relaxes 
into a pure state with all spins down. In this case, the 
symmetry between the two ordered pure states is broken 
globally. Interesting physics happens when, instead of 
a global external bias, the symmetry between the pure 
states is broken locally by the dynamics itself. In this 
paper, we investigate the effect of this particular kinetic 
disorder and show that this system also gives rise to a 
very slow dynamics. In particular, it gives rise to inverse 
logarithmic relaxation (ILR) of magnetization, very sim- 
ilar to the density compaction in granular systems. Our 
study therefore suggests that the ILR is a very robust 
phenomenon and not just limited to granular systems or 
specific types of kinetic disorders. 

The paper is organized as follows. In Section-II, we de- 
fine our model precisely and summarize the main results. 
In Section-Ill, we establish the connection between our 
spin model and a lattice model of granular compaction. 
We also show that our model can be viewed as a new one 
dimensional reaction-diffusion model when the dynam- 
ics is described in terms of the kinks between domains 
of opposite phases. In Section-IV, we derive some exact 
results. In section-V, we analyze the dynamics via an in- 
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dependent interval approximation (HA). In Section- VI, 
we argue that the IIA results become exact in the limit 
when the initial volume fraction of one of the phases is 
small. In Section- VII, we extend this argument to other 
volume fractions as well. In Section- VIII, we present a 
heuristic approach to support the IIA results. Finally we 
conclude in Section-IX. 

II. THE MODEL AND THE OVERVIEW OF 
RESULTS 

We consider a simple Ising spin chain with spins Si — 
±1. Starting from a given initial configuration the system 
evolves by single spin flip continuous time dynamics. The 
rate of flipping of a given spin depends on the its neigh- 
boring spins. We denote the rate of spin flip Si — > —Si 
by W(Si, Si-i, Si+i) where Si-i and Si+i are the two 
neighboring spins. In our model the rates are specified 
as follows: 



W (+■++) = ]¥(-■—) 

W {+■-+) = W {+■+-) 

W(+;—) = l 
W(-;++) = a 



= 

= W(-,+-)=W(-,-+) = ~ 

(1) 



Note that the case a = 1 corresponds to the usual zero 
temperature Glauber dynamics || which preserves the 
symmetry between the up and down phases. However, if 
a < 1, the flipping of a down spin sandwiched between 
two up spins is not as likely as the nipping of an up spin 
sandwiched between two down spins. Thus a < 1 clearly 
breaks the symmetry between the up and down phases. 
However this symmetry is broken only dynamically, i.e. 
not everywhere but only at at the location of the triplets 
(+,—,+). Thus the isolated '— ' spins (surrounded on 
both sides by a '+') tend to block the coalescence of '+' 
domains and locally favor the '— ' spins. One can argue 
that the asymptotic dynamics at late times is similar for 
any a < 1. In other words, a = is an attractive fixed 
point. We therefore restrict ourselves only to the case 
a = 0. 

To see the effect of the local dynamical constraint more 
precisely, we derive (following Glauber's calculation for 
a = 1 the exact evolution equation for the magneti- 
zation m(t) = (Si) for the a = case: 

^<Si) = -2P(l,-l,l), (2) 

where P((Ji_i, Ci, cr i+ i)(t) = ((l+<7 % -iS l -i)(l+a l S i )(l + 
CTi+iiSj+i)) /8 is the three point probability to find the se- 
quence of spins (<Xj_i, 0j, 0"j_|_i) about the site i and we 
have used translational invariance. Hence P(l, —1, l)(i) 
denotes the probability of the occurrence of the triplet 
'H h' at time t. We note that for the case a = 1, 



d(Si)/dt — 0, indicating that the magnetization does 
not evolve with time. In our case, due to the triplet 

defects 'H h', the average magnetization decays with 

time. If L±(t) denote the fractions of '+' and ' — ' spins, 
then using L± = (1 ± m)/2 wc find from Eq. (0), 
dL±/dt — ^fRi(t) where Ri(t) is the number density of 

the triplets of type '(H h)' per unit length, clearly 

exhibiting the asymmetry generated by the '(H (-)' 

triplets. We also note that unlike the a — 1 case, the 
evolution equation (j^) for the single point correlation 
function involves two and three point correlations (via 
i?i(t)). This hierarchy makes an exact solution difficult 
for a = 0. 

It is useful at this point to summarize the main results 
obtained in this paper. Let us first highlight the contrast 
between the a = 1 (no kinetic disorder) and the a = 
(with kinetic disorder) cases. For a = 1, due to the pre- 
served symmetry between the up and down phases at all 
times, the average domain size of both '+' and ' — ' do- 
mains grow as l±(t) ~ t 1 / 2 at late times 0. Thus the 
average magnetization m(t) = (1+ — Z_)/(Z+ + Z_) is a 
constant of motion ^]) and stays fixed at its initial value. 
In particular, if we start from an initial state where the 
spins are random (infinite temperature), the initial aver- 
age magnetization is zero and stays at all subsequent 
times. In contrast, for a = case where the symmetry is 
dynamically broken, we find the following results. 

(1) For a = 0, while domains of both phases continue 
to grow with time, they have different growth laws. The 
average domain sizes of the + phases (denoted by l+(t)) 
and — phases (denoted by l-(t)) coarsen at late times in 
the following manner: (i) l+(t) ~ \fwt at late times and 
(ii) l-(t) ~ t 1 ' 2 log(W) where & is a number depending on 
the initial volume fraction of the '+' spins which we will 
calculate explicitly (see below). In fact, the main result 
we show below is that the ratio of the two length scales 
behaves at late times as, 

I- _ 1 log(&*Ao) 1 

Z7 ~ c log b 1 {6) 

where e is the initial volume fraction (need not be small) 
and to is some initial time after which scaling starts hold- 
ing. The equation (|J) explicitly reflects the effect of bro- 
ken symmetry. Thus due to the dynamically generated 
local bias, the '— ' domains grow slightly faster than the 
'+' domains. We also point out that in contrast to the 
spin models studied in the context of glassy systems , 
the dynamics in our model does not freeze at zero tem- 
perature, rather the domains coarsen indefinitely in an 
infinite system. 

(2) Consequently, the magnetization m(t) = (l + — 

+ 1-) decays at late times as, 

, \ . 2elog b 
m(t) = -l + - — „ 6 , . (4 
v ' log(W/to) 1 ' 
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(3) Evidently the number of '+' domains is same as 
the number of of '— ' domains since they alternate on the 
1-d lattice. The density of domains of either '+' or '— ' 
type per unit length: N(t) =!/[£_ + 1+] decays as, 



N(t ) 



to log ft 
t log(W)' 



(5) 



We also compute explicitly two other densities that 
play a somewhat central role in our analysis: 

(4) n(t): Given that a '— ' domain has occurred, the 
probability that it is of length 1. We show below that 



n(t) 



V~tlog(bt/t ) 



(6) 



Note that the amplitude ^/tt is universal and independent 
of volume fraction. 

(5) pi{t): Given that a '+' domain has occurred, the 
probability that it is of length 1. We show that 



Pl(t) 



1 



1 



2t t\og(bt/t ) 



(7) 



Once again the amplitude 1 of the correction term is also 
universal and independent of volume fraction. 

The constant b in the above equations can also be com- 
puted exactly, and we show that, 



b — exp 



ri (t )Vto 



(8) 



III. CONNECTION TO GRANULAR 
COMPACTION AND REACTION DIFFUSION 
SYSTEMS 

We now establish a one to one mapping between our 
spin model (with a = 0) and a simple lattice model 
of granular compaction. Let us consider a (1 + 1)- 
dimensional granular packing where the grains are repre- 
sented by square blocks. The pack consists of horizontal 
layers consisting of blocks and voids (see Fig 1). We fo- 
cus on the 'active' layer, i.e., the first horizontal layer 
that is not fully compact with blocks as we go up from 
the bottom of the pile. Below this active layer, all layers 
are compact and remain compact under vertical tapping, 
i.e., their dynamics is completely frozen. All the activi- 
ties take place in or above the 'active' layer. We identify 
this 'active' layer with the one dimensional lattice of the 
spin model. This active layer consists of sequences of 
blocks (particles) and voids (holes). We identify a unit 
block or a particle in the active layer as a '— ' spin of our 
spin model. Similarly a hole is identified as a '+' spin of 
the spin model. 

As the system is tapped vertically, the particles in the 
active layer can undergo the following moves: (1) In the 



interior of a row of consecutive particles in the active 
layer, there is no effect of tapping as the system is com- 
pletely jammed there. The effect of vertical tapping is 
felt only at the edges of a particle cluster. The parti- 
cles at the edges, under tapping, can move up or 'tap' 
up from the active layer to the layers above the active 
layer (see Fig. 1). However if the cluster consists of an 
isolated particle sandwiched between two holes, it has no 
place to move up, so it stays at its original location. (2) 
Tapping also can 'roll' off a particle residing in a layer 
above the active layer, into the active layer, at the edges 
of the supporting cluster (see Fig. 1). 

Let us now see what the rates of spin flips in Eq. (|l|) in 
the spin model imply for the compaction model. The rate 
VF(+; ++) = implies that if we have three consecutive 
holes in the active layer, then a particle can not be de- 
posited in the middle site under tapping. This is because 
a new particle can appear into the active layer only at the 
edges of the clusters of particles, but not in the middle of 

a cluster of voids. The rate W(— ; ) = implies that 

if we have three consecutive particles in the active layer, 
the middle particle can not move up under tapping since 
it is completely jammed. The rate W(+; — h) = 1/2 sig- 
nifies that if we have three consecutive sites in the active 
layer consisting respectively of a particle, hole and a hole, 
then a particle can appear in the middle hole with rate 
1/2. This is 'rolling' off move to the right at the edge of 
a cluster, as discussed in the previous paragraph. Simi- 
larly W(+; +—) = 1/2 signifies the 'rolling' off move to 
the left. The rate W(— ; H — ) = 1/2 implies that when we 
have three consecutive sites in the active layer consisting 
respectively of a hole, particle and a particle, the mid- 
dle particle can disappear (into the layers above) with 
rate 1/2. This is the 'tapping' up move to the right from 
the edge of a cluster as discussed in the previous para- 
graph. Similarly W(— ; — h) = 1/2 signifies the rate of 
'tapping' up to the left from the edge of a cluster. The 

rate VF(+; ) = 1 indicates that if we have a sequence 

of particle, hole, particle in the active layer, the middle 
hole can be filled up with a particle with rate 1. This 
is simply the addition of two 'rolling' off rates from the 
left and the right of the middle hole. Finally the rate 
W(— ; ++) = a = implies that if we have a sequence of 
hole, particle, hole in the active layer, the middle particle 
can not disappear. This is because an isolated particle 
has no place to move or 'tap' up as discussed in the previ- 
ous paragraph. This last rate indeed breaks the particle 
hole symmetry. 

Note that the rates W(+; -+) = VF(+;H — ) = 
W{— ; H — ) = W{— ; — h) = \ correspond to the diffusion 
of domain walls in the spin model. In the compaction 
model, these are indeed the moves induced by the me- 
chanical tapping. If these rates were zero, i.e., no tap- 
ping, then the dynamics in the active layer would freeze 
after all the isolated holes are filled up and the system 
will be stuck in a metastable configuration and hence the 
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compaction will stop. It is these diffusion moves that 
lift the system out of a metastable configuration and the 
system continues to compactify, though extremely slowly. 
Identifying the '-' (+) spins with particles (holes) in the 
granular model, it is easy to see that the average mag- 
netization m in the spin model is related to the average 
particle density p(t) in the granular model via the simple 
relation, m(t) = 1 — 2p(t). Thus the result in Eq. (j|) for 
the magnetization indicates that the density will grow to 
its fully compact value 1 as 



p{t) = 1 - 



dog b 
log{bt/t ) • 



(9) 



at very late times, as observed in experiments on me- 
chanically tapped granular media Q. 

Other one dimensional lattice models, notably the 'car 
parking' model, has been previously used to explain the 
logarithmic relaxation in granular materials |§. The lo- 
cal rules for the dynamics of particles in the car park- 
ing model are, however, quite different from those in our 
model. In the parking model, the rules for the particle 
motion arc: (i) a new particle can be absorbed only at 
sites containing a hole at a rate proportional to the num- 
ber of neighboring particles and (ii) a particle evaporates 
with a small rate (thus leaving behind a hole) if it has 
one and only one neighboring hole. While this model 
also exhibits an inverse logarithmic growth of the den- 
sity, the dynamical moves (in particular the fact that the 
desorption rate is infinitesimally small) are chosen in a 
somewhat ad hoc manner. In contrast, as explained via 
the mapping detailed above, our model incorporates the 
basic minimal microscopic moves of the particles that are 
observed in the compaction process. 

In terms of the motion of the domain walls between 
'+' and '— ' phases, our model can also be considered 
a new model of two species reaction-diffusion in one di- 
mension. We note that in the case a = 1, the domain 
walls diffuse and annihilate upon contact. This corre- 
sponds to the process A + A — > @]. In the case a = 0, 
we need to distinguish between the two types of domain 
walls — h = A and H — = B. Note that by definition 
(originating from a spin configuration) the A's and P's 
always occur alternately. Here both ^4's and P's diffuse 
as before; however when an A and a B meet, they an- 
nihilate only if A is to the left of B, otherwise there is 
hard core repulsion between them. We will show below 
that this hard core repulsion between the particles is a 
relevant interaction that changes the late time dynamics 
considerably. Recently the relevance of hard core repul- 
sion between particles in reaction-diffusion systems have 
been explored in a number of contexts JlCj ] . 



IV. SOME EXACT RESULTS 

To start with we write down two exact relations which 
are derived directly from the Glauber dynamics (we shall 
see later that these exact relations are in fact respected 
by the IIA approximation). 

Let N(t) be the number of domains of either '+' or ' — ' 
spins per unit length. The density of kinks is therefore 
2N(t). Let Pi be the density of the triplets ' — I — ' per 
unit length, then clearly pi — P\/N . We note also that 
the only way in which a domain can be destroyed in an 
infinitesimal time step is by flipping an isolated '+' spin 
in a triplet ' — I — '. This gives 



dN 

~dT - ~ Pl ' 



(10) 



If P(Si) denotes the probability that the spin at site i 
takes the value Si, the evolution of P(Si) depends only 
on the rates in Eq. (0) and the three point probability 
distribution P(S'j_ 1 , Si, Si+i)- The evolution of P{Si) is 
then given by 

dP(Si) 



dt 



W(~Si, Si— i, Si+i)P(Si-i , —Si, Si+i) 



E 

S i+ i,Si- 



W(S l ;S l - 1 ,S l+1 )P(S l - 1 ,S l ,S l+1 ) (11) 



Substituting Si = 1 in the above equation we find 
dP(l) 



dt 



P(-l, -1, 1) - P(-l, 1, -1) - P(l, 1, -1) 



(12) 



where we have used the evident left to right symme- 
tries P(-l,-l,l) = P(l,-1,-1) and P(l,l,-1) = 
P(-l,l,l). We now observe that P(-l,-l,l) + 
P(l, -1, 1) = P(-l, 1) and P(-l, 1, -1) + P(l, 1, -1) = 
P(l, -1) = P(-l, 1) and thus find 

dP{\) 



dt 



= -P(l,-l,l). 



(13) 



We note that P(l,— 1,1) is simply the probability that 
at a given site the spin is a '— ' spin and its two neighbors 

are '+' spins, that is to say that their is a 'H h' defect 

at the site considered. If we now sum this equation over 
each site on an interval of unit length on the lattice and 

recall that Pi is the density of the 'H K triplets per 

unit length, we obtain 

where L + is the fraction of the '+' spins. Note that 
L + + L- = 1, where L_ is the fraction of the '— ' spins. 
From the relation Eq. (|l4|) we obtain Eq. (||) for the 
evolution of the average magnetization m. Physically it 
is easy to see the origin of Eq. ( |l4| ) as, on an average, 
the fraction of '+' spins can decrease only due to the 
blockage by l H h' triplets. 
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Writing Eq. (|14|) in terms of r% = R\ /N and the aver- 
age length of the '+' domains l + = L + /N we obtain 



N 
N 



l+ + ri 



(15) 



where x = dx/dt. This equation can be integrated start- 
ing at some arbitrary time to to give 



N(t) l + (t Q ) 



N(t ) l+(t) 



exp 



ri(t') 

i+(t'y 



dt' 



(16) 



Furthermore, if the volume fraction of the '+' phase is, 
L + (t ) = e, then, using the relation, N(t) — 1/[Z_ (t) + 
l+{t)} in Eq. (pi), we find, 



MO 
hit) 



1 



■ exp 



l+(t'Y 



dt' 



(17) 



clearly showing that the ratio l^(t)/l+(t) is growing due 

to the presence of the triplets 'H K. Note that the 

asymmetry between the growth of '— ' and '+' domains 

is evident due to the presence of the triplet defects 'H h 1 

with density R\ — r\N . 

All the results presented above are exact. To derive 
the late time behavior of the model we first consider be- 
low the IIA. We solve the IIA equations self-consistently 
and show that the IIA precisely predicts the results men- 
tioned in Section II. Besides we shall argue that in the 
case where the initial volume fraction of the '+' domains 
is small, i.e., e < 1, correlations do not develop between 
the domains if no correlations are present in the initial 
conditions and hence the IIA is exact to leading order in 
e. At the end we present a very simple heuristic argument 
which is also in agreement with these results. 



V. IIA ANALYSIS 

In this section we consider the IIA where correlations 
between neighboring domains are neglected. The IIA was 
used previously for the a — 1 case jn]] yielding results 
in agreement, qualitatively as well as quantitatively to a 
fair degree of accuracy, with the exact results available 
HQ- Let P n (t) and R n (t) denote respectively the num- 
ber density of '+' and '— ' domains of length n at time 

t. Note that Ri(t) is the density of the triplet 'H h ! 

as before. Let N(t) — J2 n ^ n ~ S n ^« denote the do- 
main density of '+' or '— ' spins. Also the fractions L±(t) 
of '+' and '— ' spins are given by, L + (t) — ^2 n nP n and 
L -{t) = Y,„ nR n with L+(t) + L-(t) = 1. During an 
infinitesimal time step At, P n (t) evolves as: 



P n (t + At) = P n (t) - AtP n (t) - AtP n (t) 



1 



Ri{t) 



AtP n+ i(t) + A«P„_i(t) 



1 - 



N(t) 



(18) 



The right hand side of the above equation includes the 
various loss and gain terms. The second and the third 
terms describe respectively the loss due to the hopping 
inward and hopping outward of the domain walls at the 
two ends of a '+' domain of size n. An outward hop can 
occur provided the neighboring domain in the direction of 
the hop is not an isolated '— ' spin and this is ensured by 
the prefactor (1 — Ri/N) in the third term. The fourth 
and the last term describe similarly the corresponding 
gains. One can similarly write down the evolution equa- 
tion for the i?„(i)'s. During an infinitesimal time step 
At, R n (t) evolves as: 



R n {t + At) = R n (t) - AtR n (t) - AtR n (t) 
- 2AtP 1 (t)^l + AtR n+1 (t) 



1 - 



Px(t) 



N(t) 



N(t) 



+ AtR n ^(t) 



Pi(t) 



N(t) 



Pi(t) 
N 2 (t) 



71-2 



Y,Ri(t)R n - t -i(t) :n>2 (19) 



and 



Rt(t + At) = R t (t) - AtR^t) 
-2AtP 1 (t)ff 4 



Pi(t) 

- im. 

AtR 2 (t). 



(20) 



The negative (loss) terms in Eq. (|19|) for domains of size 
n > 2 may be understood as follows. A domain of length 
n may be lost by the domain wall at either end jumping 
inwards with rate 1/2. This term is the second term in 
Eq. ( |l9| ) and as there are two domain walls we have a 
factor of 2. A domain of length n may also be lost by a 
domain wall hopping outwards. This happens with rate 
1/2 if the neighboring domain is not a triplet ' — I — '. 
The third term of Eq. (|n]) corresponds to this event, the 
factor (1 — Pi(t)/N(t)) is the probabilty of the absence 
of a triplet 1 — | — ' as a neighboring domain. There is 
again a factor of 2 coming from the fact that there are 
two domain walls. However if a neighboring domain is 
of type ' — I — ' the outward jump towards this domain 
occurs with rate 1 (as the central + spin flips with rate 
1). The term corresponding to these two events (from 
the right and left domain walls) is the fourth term in Eq. 
( |l9| ) . The two first gain terms come from identical argu- 
ments and the last convolution term represents domain 
coalescence, where a domain of length n is formed with 
rate 1 from two '— ' domains of length % and n — i — 1 
(where 1 < i < n — 2) if they are separated by a ' — | — ' 
triplet. Eq. ( p0[ ) is obtained in a similar fashion with 
the exception that the hard core repulsion generates a 
reflecting boundary condition. Taking the limit At — > 
in the above equations we obtain the IIA equations for 
the evolution of the domain densities 
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dP-n R\ 

—j— = P n +1 + P-n-l — 2-P n + -jr{ P-n ~ Pn-1) 



(21) 



for all n > 1 with P Q = (absorbing boundary condition) 
and 

dR n Pi 

= Rn+l + Rn-l — 2R n — -j^(Rn + Rn-l) 



' N 2 



n-2 



R%Rn-%-i\ n > 2 



dRi Pi 

— = R 2 - Rl -- Rl , 



(22) 



It is somewhat convenient to use the normalized vari- 
ables p n — P n /N and r n = R n /N . The average domain 
lengths are then given by l+(t) = ^np n and l-(t) = 
nr n and the domain density N(t) = l/(l+(t) + l-(t)). 
In terms of these normalized variables the IIA equations 
are given by, 



dt 



: Pn+l + Pn-1 - 2p n + ri (p n - Pn-l) + PlPn (23) 



for all n > 1 with pq = (absorbing boundary condition) 
and 

dv 

—jj- = r n+1 + r n -i - 2r n - p\r n _ x 



+Pi 2J nr n -i-i', n>2 



i=l 



dr\ 
~dt 



r 2 - r x . 



(24) 



It is easy to check that the normalization condition, 
T^Pn — J2 r n — 1 is satisfied by these two equations. 
We note that Eq. (2^) or more clearly its unnormalized 
version in Eq. (|2l|), i.e., P n = P n+1 + P„_ x - 2P n + 
fi{P n — Pn-i), just represents the motion of a random 
walker on the positive side of a 1-d lattice with a sink at 
the origin (Pq = 0) and a time dependent drift term (pro- 
portional to ri). To calculate N(t) using Eq. (|ilf), we 
need to evaluate two quantities from the IIA equations: 
(i) n(t) - Ri/N and (ii) /+(*) = Y. n Pn- 

The two IIA equations above are coupled nonlinear 
equations with infinite number of variables and hence 
exact solution of Eqs. (|2^) and (^4|) are difficult. Our 
approach will be a combination of a scaling assumption 
and then rechecking this assumption for self-consistency. 
Consider first the r n equation, i.e. Eq. (p4|). On the 
right hand side, we will first ignore the diffusion term, 
solve for the rest and show that indeed neglecting the 
diffusion term was justified in the first place. This is 
self-consistency. Ignoring the diffusion term, we have the 
following equation, 



dr n 
dt 



Pi 



n-2 

£< 

i=l 



-1 - r n _i 



(25) 



with the reflecting boundary condition, r 2 — r\. It is 
now easy to see that Eq. (^J) admits a scaling solu- 
tion, r n (t) rj X(t) exp[— n\(t)}, where X(t) = —pi(t)\(t). 
Using the exact relation, dN/dt — —piN, we get, 
A(t) = X(t )N(t)/N(t a ). Note that we still do not 
know what N(t) is. Now let us substitute this solu- 
tion to estimate the diffusion term that had been ne- 
glected in the first place. Clearly, the diffusion term, 
Tdiff = f*n+i +f"n-l — 2r n ~ 0(X 3 (t)), whereas the other 
terms (for example the left hand side of Eq. (|24|)) typi- 
cally scale as ~ 0(X(t)) ~ 0(pi(t)X(t)). Thus, in order 
to be self-consistent in neglecting the diffusion term, we 
need to havepi(t) > X 2 (t) — N 2 (t). We will see that this 
condition is actually satisfied once we derive the expres- 
sion for N(t). This just means that the diffusion terms 
only contribute to the corrections to the leading scaling 
behavior. 

From the above analysis, we find to leading order 
for large i, n sa X(t) = X(t )N(t)/N(t ), i.e., n(t) pa 
ri(to)N(t)/N(to). We now need to evaluate the other re- 
maining quantity, l+(t) = J~] np n . For this we now turn 
to the p n equation, Eq. (|23|). In this equation, we will 
again first ignore the drift term ri(p n — Pn-i) solve for 
the rest and check that indeed the neglect of the drift 
term was justified. Ignoring the drift term, we get, 



dp n 
dt 



—7T = Pn+l + Pn-1 - 2p n + pip n . 



(26) 



with the absorbing boundary condition, pq = 0. This 
equation can be solved exactly. Indeed it also admits a 
scaling solution, p n (t) = t~ x l 2 /(nt^ 1 / 2 ), where the scal- 
ing function (normalized to unity) is given by, f(x) — 
I exp(— x 2 /A). Now let us estimate the drift term that 
was neglected. Clearly the drift term, r\(p n — p n -i) ~ 
0(ri/t) ~ 0(N(t)/t) since r\ ~ N(t) from previous 
paragraph. The other terms in the Eq. ( |23| ) (for ex- 
ample the left hand side of Eq. (p3|)) is of order, t~ 3 / 2 
at late times. Thus, for self-consistency in neglecting 
the drift term, we need to show that £~ 3 / 2 ^> N(t)/t. 
We will again see that this condition is indeed also sat- 
isfied once we obtain the expression for N(t). From 
this form of p n (t), we thus obtain, to leading order 
for large t, l+(t) = J2 n Pn ~ t 1 ^ 2 xf(x)dx. Using 
f(x) = I exp(— x 2 /A), and doing the integral we finally 
find, l+(t) ~ \pnt for large t. 

Using these two results (i) n(t) = ri(t )A^(/j)/A r (t ) 
and (ii) l+(t) = \f%i in the exact equation Eq. (|l6|) and 
differentiating with respect to t, we find a differential 
equation for N(t), 



±(ViN(t))=- 



dt 



rijtg) 
V^N(t Q 



■N 2 (t). 



(27) 



Introducing the dimcnsionless variable 



G 



S(t) 



rj* N(t) 



t N{to) 



(28) 



in Eq. (U7V) we find 



dS S 2 

^ = - l0g(6) T 



(29) 



where log(6) = 



' , i(*o)\/*o ' 



Integrating Eq. (09]) we find 



S(t) 



log(&) 
log(W/t ) ' 



(30) 



we thus obtain the result 



N(t) 



N(to 



t ]og(bt/t y 



(31) 



Substituting this result in the expression, 7*i(t) 
ri(t )N(t)/N(t ), we get 

n(i) 



Vtlog(6i/i )" 



(32) 



Next we use the late time result Eq. ( |3l| ) in the exact 
relation Eq. (fLOh and find, 



(33) 



£>i = . 

^ 2i tlog(bt/t ) 



Let us check the two self-consistency conditions, (a) 
pi(t) > A 2 (i) - N 2 (t) and (b) t~ 3 / 2 > N(t)/t. Using 
the expression for pi from Eq. ( |33| ) and that of N(t) 
from Eq. (|3l|), it is immediately evident that indeed 
these two conditions are satisfied for large t. Thus our 
whole approach has been completely self-consistent and 
the IIA results are precisely those mentioned in the ab- 
stract. Note also that to must be sufficiently large such 
that both scaling laws, ri(t) ~ N(t) and l+(t) ~ y/wt 
start holding for t > to. 



VI. ZERO VOLUME FRACTION LIMIT 

In this section, we show that the IIA results essen- 
tially become exact in the zero volume fraction limit of 
the '+' phase, i.e., in the limit e — > 0. Suppose we start 
from an initial condition such that, r„(0) = e(l — e) n and 
Pn(O) = 8 n ,i- This means that in the initial condition, the 
average length of the ' — ' domains, Z_(0) ~ 1/e, whereas 
Z+(0) = 1. Thus the '— ' domains are typically much 
larger than the '+' domains, in the limit e — * 0. Also 
initially all the domains are completely uncorrelated. So 
the picture is as follows. We have little droplets of '+' 
phase in a sea of '— ' phase. Besides, one can also com- 
pute the initial density of domains of either '+' or '— ' 
types. It is given by, N(0) = l/e(l — e) ~ 1/e to leading 
order in e. 



Now let us consider the time evolution of the system 
starting from this initial condition. As time increases, 
the '+' domains will certainly grow in size. But a typical 
'+' domain will disappear (via the absorbing boundary 
condition) much before encountering other '+' domains, 
i.e., before feeling the presence of the constraint due to 

triplets 'H h'. The probability of such an event is of 

order 0(e). Thus effectively, the dynamics of the system 
will proceed via eating up of the '+' domains. Hence, 
if there is no correlation between domains in the initial 
condition, the dynamics is not going to generate correla- 
tions between them. This is precisely what happens in 
the zero temperature dynamics of the q state Potts model 
in 1-d in the limit q -> 1+ US 

Thus in this limit, the evolution of the ' — ' domains is 
governed by the exact equation, 



dr n 
dt 



Pi 



n-2 



(34) 



which is same as the IIA equation Eq. without 
the diffusion term. Starting from the initial condition, 
r n (0) = e(l — e) n , one can solve the above equation for 
any t exactly to leading order in e. It turns out that, to 
leading order in e, Eq. (34) admits a solution, r n (t) = 
H(t)[l - mC*)]"" 1 ) where fi(t) = eexp[- Pl {t')dt'}. Us- 
ing once again, the exact equation dN/dt = —piN, we 
find, fi(t) = eN(t)/N(0) + 0(e 2 ) = N{t) + 0(e 2 ) where 
we have used N(0) = 1/e + 0(e 2 ). Thus we get, 



ri(t) =m(*) = N(t) + 0(e 2 ). 



(35) 



Now, let us consider the evolution of the '+' domains. 
Since a typical '+' domain never encounters (to leading 
order in e) any other '+' domain and hence does not feel 
the constraint due to ri's, the effective dynamics of a '+' 
domain is that of a single '+' domain immersed in a sea 
of '— ' phase. Let P n denote the probability that such a 
domain is of length n. Then, to leading order in e, P„'s 
clearly evolve by the simple diffusion equation, 



dP n 

dt 



Pn+l + Pn—1 — 2P„, 



(36) 



with the absorbing boundary condition, Pq = 0. The 
normalized conditional probability, p n = P n /N with 
N = Pn, then satisfies the equation, 



d Pn 

—- = p n+ i + p n -x - 2 Pn + P i Pn 

dt 



(37) 



same as Eq. (p6|). This equation has to be solved with 
the initial condition, p n (0) = 1. It can be solved exactly. 
Without writing the explicit solution, we just mention 
the result for np n . We find that for large t and lead- 
ing order in e, 



*+(*) = £ 



n Pn 



7rt + 0(e). 



(38) 
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Using the results from Eqs. (|3|) and (|38|), i.e., (i) 
ri(t) = N(t) + 0(e 2 ) and (ii) l+{t) « V^rf + 0(e) in the 
exact equation Eq. ([!(]), we once again recover all the 
IIA results of the previous section , with b — vW 6 - 

Hence IIA becomes exact in the e — > limit. This 
is not surprising as the dynamics in this limit does not 
generate correlations if there are none in the initial con- 
dition. 

VII. OTHER VOLUME FRACTIONS 



are now shrunk to a single point when viewed from the 
length scale of the '— ' domains. The rate of occurrence 
of these points per unit length X(t) is clearly propor- 
tional to N(t), i.e., the kink density. Note that r n is 
simply the conditional probability: Given that a '— ' do- 
main has occurred, what is the probability that it is of 
length n. Now if one assumes that these punctual '+' 
domains (or the points) are distributed randomly, one 
finds that r n is simply given by the geometric distribu- 
tion, r„ = A(i)(l - A(i))™" 1 where A(i) = cN(t) for some 
constant c. We therefore find 



For finite initial volume fraction of the '+' phase, the 
IIA can not be exact since the diffusion of kinks corre- 
late the domains as time progresses, even if the domains 
had no correlations to start with. However, the volume 
fraction of the '+' phase decreases monotonically with 
time according to the exact equation Eq. (|l4|). Thus at 
very late times when the volume fraction L+(t) is very 
small, the effective fixed point picture of the system is 
very similar to the e — > limit picture, i.e., small '+' 
domains immersed in the sea of '— ' domains. The only 
difference is that the big '— ' domains may now be cor- 
related. However, it is very likely (though we can not 
prove this rigorously) that the correlations between do- 
mains are very small at very late times and therefore the 
IIA results ((l)-(4) in the abstract) become asymptoti- 
cally exact. The numerical results reported in [pi appear 
to confirm this fact. 

Actually if this late time fixed point picture is correct, 
then one can derive all the results from a very simple 
heuristic argument, presented below. 



VIII. HEURISTIC APPROACH 

The heuristic picture is as follows. Due to the facts 
that '+' domains can grow only by diffusion and the '— ' 
domains grow by diffusion and coalescence, one expects 
that at late times the ratio i+(i)/7_(t) — > 0; this is also 



n(t) wdV(t), 



(40) 



clear from Eq. (14). Therefore at late times one expects 
to find '+' domains sandwiched between much larger '— ' 
domains. Consequently the '+' domains, to a first ap- 
proximation, never encounter the 'H (-' triplets at late 

times and hence diffuse freely and annihilate via the dis- 
appearance of ' — I — ' triplets. Therefore the length of 
a +' domain behaves as a diffusion with a killing bound- 
ary condition at the origin. Solving the discrete diffusion 
equation corresponding to this picture (see Eq. (|37|)), 
one finds that at late times 



l + (t) 



nt 



(39) 



If one now views the system at the length scale of the 
'— ' domains l-(t) 3> one sees long stretches of '— ' 

domains occasionally interrupted by '+' domains which 



for late times. Furthermore, if we denote by to a large 
time after which this picture becomes valid we may write 



ri(t) w T-i (t ) 



N(t) 
N(t ) ' 



(41) 



These two expressions for r\ and l+(t) once again are 
same as obtained by a more careful analysis of the IIA 
equations and when substituted in Eq. ([l6|), they give 
the same IIA results once again. Thus the basic assump- 
tion of this heuristic picture is the '+' domains occur 
randomly, which seems like an accurate description at 
late times. 



IX. CONCLUSIONS 

In this paper we have presented detailed analytical 
studies on a simple one dimensional kinetically con- 
strained Ising model which was introduced in Ref. ||. 
The kinetic constraints in this model are local and dy- 
namically generated. The effect of these constraints was 
shown to slow down the dynamics rather drammatically. 
We have shown that the average magnetization in this 
model decays extremely slowly with time in an inverse 
logarithmic fashion to its final saturation value. This 
kind of inverse logarithmic law was observed in the be- 
havior of the density of granular material in experiments 
on granular compaction || and was also seen in numer- 
ical simulations of various lattice based and 'tetris' like 
models fl4"| . There have been some theoretical arguments 
proposing various mechanisms responsible for this slow 
compaction 15,^]. These include the free volume argu- 
ment fll5| and the argument based on an analogy with 
car parking models with the somewhat ad hoc assump- 
tion that the cars 'depark' from a lane at an infinitesimal 
rate [Bj. In contrast, in this paper we have mapped our 
kinetic Ising model directly to a lattice model of granu- 
lar compaction which incorporates the basic minimal mi- 
croscopic moves in the compaction process. The average 
magnetization m(t) in the Ising model, via this mapping, 
gets related to the density of compaction p(t) in the gran- 
ular model as p(t) = (1 — m)/2. Hence, besides having 
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nontrivial behavior and yet analytically solvable, our toy 
model of granular compaction correctly reproduces the 
inverse logarithmic time dependence seen in the experi- 
ments ||] and thereby proposes a new and entirely differ- 
ent mechanism for this slow compaction, quite different 
from the previous models such as the car parking models. 
It is also interesting to note a study of compaction in the 
Tetris model [[t6| shows that at late times the activity of 
the system, leading to compaction, occurs at boundaries 
between domains which can be identified in the system. 
The image of compaction as a kinetically hindered coars- 
ening process thus appears to be quite robust. 

From a somewhat broader perspective, our work ad- 
dresses a general question: what is the effect of kinet- 
ically generated disorders on the coarsening dynamics 
in domain growth problems? In the present work we 
have studied a specific type of kinetic disorder, namely 
a dynamically generated local magnetic field. This field 
acts locally on the topological defects responsible for the 
coarsening process (in this case simple domain walls). 
Our study suggests that such kinetic disorders, while 
slowing down the dynamics drastically, do not altogether 
inhibit the coarsening process as found in other con- 
strained kinetic Ising models ||. The domain growth 
problems are rather common and occur in various physi- 
cal systems (|. Our work, therefore, opens up the possi- 
bility of studying the slowing down in coarsening dynam- 
ics due to kinetic disorders in many of these systems. For 
example, it would be interesting to study the effect of dy- 
namically generated local fields in higher dimensionional 
Ising models, in 0{n) vector spin models and in liquid 
crystals, to mention a few. It is possible to have other 
types of local kinetic disorders than the one studied here 
and it would also be interesting to study their effect in 
coarsening systems. 

We thank P. Grassberger for his earlier collaboration 
in this study, and also M. Barma for useful discussions. 
We acknowledge interesting exchanges with G. Odor, H. 
Hinrichsen and G. Schiitz. 
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Figure Captions 

Fig. 1. Picture of rolling off (left domain) and tapping up (right domain) in the granular interpretation of the spin 
model, the solid squares represent particles. 
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